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Mutual Impedances of Grounded Circuits 

By GEORGE A. CAMPBELL 

Synopsis: Formulas are derived for the direct -current mutual resistance 
and inductance between circuits grounded at the surface of the earth. For 
circuits composed of straight filaments, the mutual inductance is reduced 
to known Neumann integrals which involve only comparatively simple 
expressions for the case of horizontal, coplanar conductors above, below 
or on the surface of the earth. Numerical values for these integrals may 
be readily obtained from new and accurate graphs for straight filaments 
which meet at a point or start from a common perpendicular. It is shown 
that these new results supply a useful first approximation to the actual 
alternating-current mutual impedance of grounded circuits, when the fre- 
quency and extent of the circuits are not larger than occur in many prac- 
tical applications. 

1. Introduction 

THE important discovery of the possibility of using the earth 
as the return conductor for electric telegraphic communication 
was announced by Steinheil in the Comptes Renins of September 10, 
1838, and throughout the entire development of telegraphy grounded 
circuits have been extensively employed. Considering the extensive 
application of such a capital discovery extending over a period of 
85 years, it is surprising that so little is known quantitatively about 
grounded circuits. We have, however, long known that conditions 
are not of the extreme simplicity pictured under the early view that 
the earth acts as a reservoir presenting no resistance to the return 
current and introducing no interference between parallel returns. 
This view was expressed in 1857 by Bakewell as follows: "There is 
no mingling of currents, the electric current of each battery being 
kept as distinct as if separate wires were used both for the trans- 
mitted and the return currents. It would indeed be as impossible 
for the separate currents transmitted from the two batteries to be 
mingled together as it would be for the written contents of two letters 
enclosed in the same mail bag to intermix." 

Measurements made a few years ago of the mutual impedances 
between grounded circuits which are restricted to a territory six miles 
square, at frequencies of 25 to 60 cycles per second, showed that within 
10 per cent, the mutual reactance increased in the same ratio as the 
frequency. It was inferred that the effective inductance under the 
Conditions of these tests was approximately the same as for direct 
current, or in other words, the incomplete penetration of the alter- 
nating currents into, the earth was not of controlling importance in 
tests upon this scale. 
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This led to my making a theoretical investigation of the mutual 
inductance between direct-current grounded circuits which did, in 
fact, show that the calculated numerical results are in reasonable 
agreement with these actual experimental data. It is the purpose of 
this paper to describe this work; the mathematical discussion of the 
theoretical .corrections for the incomplete penetration of alternating 
currents into the earth will form the subject of another paper. 

2. Distribution of Current, Potential and Magnetic Force 
with Direct-Current Earth Return Flow 

On the assumption of an infinite earth of uniform resistivity the 
lines of flow and the equipotential surfaces for a direct current / 
entering the earth at a point source at A and leaving the earth at a 
point sink at B, both A and B being on the surface of the earth, 
assumed flat, and the distance AB = 2b, are given by the equations, 1 
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where C is the total current flowing in the earth, from the source at A 
to the sink at B, outside the current sheet of revolution defined by (1) ; 
and V is the potential, with respect to the midplane, upon the equi- 
potential surface of revolution defined by (2). 

These equipotential lines and stream lines are identical with the 
equipotential lines and lines of force for a uniformly magnetized 
filament. The formulas may be checked by regarding the return 
flow from A to B as being due to the superposition of two flows, a 

»The coordinates used in this paper (*i, y, a), (xs, y..«0. (x, y, z) and (r u lf 0), 
(f2 02, <t>), (r, 0, <i>) are rectangular and spherical coordinates with origins at A. B 
and the midpoint of AB, respectively, the direction AB being the polar axis or posi- 
tive x-axis in all cases, z being vertical, and 4> being measured from the eafth'i 
surface in thejplane perpendicular to AB. 
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return flow of direct current / from the point A to some infinitely 
distant point and a second return flow of direct current / from this 
infinitely distant point to the point B. For these component flows 
the current diverges from A or converges towards B radially and 
with equal intensity in all directions in the earth; the total current 
for one of the component flows flowing through any surface in the 
earth will thus be equal to 1/2* times the solid angle subtended at 
A or B, respectively, by the boundary of the surface, s : nce the entire 
solid angle filled by the earth at a point on the surface is 2tt. The 
total radial flow from A through the lower half of the circular cone 
having its axis in AB, the elements of the cone making the angle 
0i with AB, is %I(l— cos 0i); similarly, the total radial flow toward 
B through the lower half of a cone with the angle ir — 2 will be 
i/(l-f-cos 2 ). For the combined superposed flows the total current 
flowing through the semicircle in which the cones intersect is the 
sum of these two values or \I{2 — cos 0i+cos 2 ), from which (1) 
is immediately obtained, since the total current flowing in the earth 
from A to B is /. This assumes that the semicircle lies between A 
and B, but the same formula holds for the entire current sheet of 
revolution. The lines of flow in the earth are symmetrical about AB 
and lie in planes through AB, since, in the earth, both component 
flows are symmetrical about AB. 

For the component flows the equipotential surfaces are hemispherical 
and, since the resistance of a hemispherical shell of radius r, thickness 
dr, is pdr/2-irr 2 , the potentials at distances r x or r 2 from ^4 or B, re- 
ferred to the potential at infinity, are Ip/2irry or — Ip/2irr 2 , respec- 
tively, from which equation (2) follows by addition. 

Fig. 1 accurately reproduces the flow and equipotential lines as 
given by formulas (1) and (2). At the midpoint of a line of flow its 
distance from each electrode is ri = r 2 = &//C and it may be shown 
that every other point of a line of flow is at a still shorter distance 
from the nearer electrode. It follows, for example, that less than 
1/10 of the total current reaches, in its flow through the earth, any 
point lying at a distance greater than 5AB from the line AB connect- 
ing the electrodes. 

If a uniform radial flow of current / in the horizon plane con- 
verging on the point A is combined with the uniform radial flow 
in the earth outward from A , we have a closed flow which is sym- 
metrical about the vertical axis through A. Below the horizon 
plane the magnetic lines of force will be horizontal circles and the 
magnetic force at any point distant r x from A, a being the angle 
included between r x and the nadir, is H = 2I(l — cos a)/(r isin a) 
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= 21 fr 1 tan (a/2). Above the horizon plane there is no magnetic 
field, since any magnetic lines of force are, by symmetry, horizontal 
circles and the intensity is zero, since there is no current threading 
any horizontal circle above the surface of the earth. 

Superposing this closed flow and a similar closed flow through B 
from the earth to the horizon plane, we obtain a closed flow from A 




io'-° 



Fig. 1— Flow and cquipotential lines on the earth's surface for an earth return flow 
from A to B. C/I is the fraction of the current flowing in the earth outside the 
flow surface of revolution; V/I is the resistance to the flow of the portion of the 
earth lying between the equipotential surface of revolution and the mid or zero po- 
tential plane, if the earth's resistivity is p = 2ir. The flow and equipotential lines 
through each point on the dotted curve are perpendicular and parallel, respectively, 

to AB 

to B in the earth and back from B to A in the horizon plane. The 
magnetic field for this closed flow, being the sum of the magnetic 
fields for the component flows, will be zero above the horizon plane, 
while below it will consist of lines of force in horizontal planes. This 
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result also applies to any closed flow which does not extend above 
the horizon plane and may be resolved into any number of com- 
ponent flows, each of which is radially symmetrical about a vertical 
axis. 

3. Mutual Resistance of Grounded Circuits 

By definition, if e.m.f.'s E and e in grounded conductors AB and 
ab produce the currents / and i = in the conductors, the mutual 
impedance between the two conductors is e/I. In the present case 
we are dealing with direct current and thus the mutual impedance 
is a mutual resistance Q, and by (2) its value is ■ 



Ab 



Ba ^ Bb) 



d? (\ 



^ 2tt \Aa 

= 2rJ J dSds \R, 

P r f - 2d Udu +dVdv+d Wdw 

2ttJ 



dSds 



R" 



cos (0 X — e) cos (0 2 — «) 



ds. 



(3) 



The third form of (3) shows that the mutual resistance falls off 
as the inverse third power of the distance between grounded circuits 
when this distance has become large compared with the length of 
these circuits between grounding points. 

The first form of (3) shows that the mutual resistance between 
grounded circuits does not depend upon the location of the conductors 
but only upon the location of the terminal grounding points A, B, a, b. 

The mutual resistance for the case p = 2ir is obtained from Fig. 1 
by taking the value of V/I at the point corresponding to a reduced 
by its value at the point corresponding to b; if b is anywhere on the 
center line, for which V/I = 0, the diagram gives directly the value of 
the mutual resistance. Employing ordinary units the diagram gives 
the mutual resistance directly in ohms if AB = l mile and the earth 
has a resistivity of about one million ohms per centimeter cube (more 
exactly 1.011X10 6 ) which is its actual order of magnitude. 

2 In addition to the earlier notation there are employed in the different expressions 
for formula (3), and also in formula (5) below, the following: R is the distance be- 
tween two elements dS and ds of any two paths extending from A to B and a to b; 
the rectangular projections of these elements along and perpendicular to R are dU, 
dV, dM^and du, dv, dw, the two sets being parallel and with the same positive direc- 
tions; 0i, 02, t are the angles which r u ri and ds make with AB, when the path ab 
lies in a plane with A and B. 
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4. Neumann Integrals for Return Flows 

The required mutual inductances of grounded circuits will be found 
by means of the Neumann integral 

N = ff f ("(cos e/r)dl dS di ds 

extended over every current filament in both flows. Since the earth 
return portions of the two flows are independent of the flows in the 
arbitrarily located conductors on the earth's surface, it is convenient 
to divide the Neumann integral into four partial integrals which in- 
volve either no return flow, one return flow or both return flows ac- 
cording to the following formula 3 

N(jr- S) U- ») = NjTx— Njt. - Nsjt-\- Ns* 

= Njr Jr -(h+h-l)A, by Table I, 

= Njf Je . (4) 

Checking the entries of Table I may be accomplished without per- 
forming more than two integrations. It will be convenient to make 
the integrals somewhat more general than is required in checking 
the table and find Njr s and Njt'a where J? is any flow in space from 
■A to B, which need not be coplanar points with the terminals a and b 
of •», and -X" is any flow in a plane parallel to the horizon plane be- 
tween terminal points A' and B' . 

Consider first the part of a space return flow * which is radial from 
a in connection with an element dS on any filament of current di 
of a flow •& from A to B. The component dx of dS along the line x 
from a to ds is the only component which need be considered, since by 
symmetry the normal component contributes nothing to the Neumann 
integral. As the total radial flow is to be taken equal to unity, the 
amount flowing out through a ring, taken as the volume element, 
lying between the spheres of radii 5 and s-\-ds and between the cones 
making angles 6 and 6-\-d8 with x will be \ sin 6 dd. If this ring lies 
at a distance r from dS the Neumann integral will be 

N= dx I ds I » ,r* = x i + s* — 2xs cos d, 

\J Aa xVo JV|*_,| 

3 Each term indicates the Neumann integral for the pair of flows designated by 
the script letter subscripts, as explained in the note accompanying Table I. Both 
{JT—G) and (x—e) are arbitrary flows on the earth's surface closed by earth return 
flows from A to B and from a to b, respectively. 



MUTUAL IMPEDANCES OF GROUNDED CIRCUITS 



i r Ba 
= U dx = \{Ba-Aa). 



ds 
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Table I 

Entries are the value of k in the formula N = kA for the Neumann integral between 

the specified flows, A= -Aa+Ab+Ba-Bb, or 2/15 if A =a, B = b, and points A 

B, a, b are all on the earth's surface. 
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* Not of the assumed form and in general infinite. 

t Each type of flow is designated by a script letter. Any flow in the surface of 
the earth, assumed to be a plane, is designated by <X~, there being no restriction on 
the number of filaments of current but each filament must start from a common 
point, A, and extends to another common point, B. A special case of JTis ^"the 
horizon return flow made up of two superposed uniform radial flows, from A to 
every point of the horizon and back to B. The space return flow, S, is made up 
of two superposed uniform radial flows, one outward in all directions from the point 
A and the other inward from all directions toward the point B. The earth return 
and air return flows, £ and J^ are similar except that the flows are uniformly dis- 
tributed over all directions in the earth and air, respectively. The nadir return 
and zenith return flows, 9? and .2", consist of two infinite vertical filaments from A 
and back to B by way of the nadir and zenith, respectively. Small script letters 
indicate similar types of flow with any independent terminal points, a and b. Dif- 
ferences such as {Jf—S) designate closed flows; thus {JC—G) designates any flow 
on the earth's surface from A to B where it enters the earth and after spreading 
out uniformly through the earth returns to the terminal A, thus closing the flow. 
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To this is to be added the corresponding expression \{Ab — Bb) 
for the radial flow converging on b, giving the result £A. As the path 
of the line of flow between A and B does not enter into the result, 
it is immaterial whether the flow is confined to a single filament or 
is spread out in any way whatsoever in space, provided only all 
stream lines extend from A to B, as assumed for &. Thus 

To find Njit'a let r and s be the distances of any element dS of a 
line of flow forming a part of <X~' from a and from any element of a 
plane radial flow from a, respectively, the projections of r and s on 
either of the planes being* and y, which include the angle 0; s 2 = y 2 -\-z 2 , 
r 2 = x 2 -\-z 2 . The component of dS parallel to x will be dx and this is 
the only component which need be considered, since, on account of 
the symmetry of the radial flow, the normal component in the plane 
of flow contributes nothing to the integral. 



dx / -YT 

..j Jo Jo x z + 



x — y cos <p y dy d(f> 



l -\-y 2 — 2xy cos 2irs 

= hf Aa il ds f Q ( 1 + «»+/-2xycosJ^ 

_ J_ pBa dx f M , p . _ 1 (x 2 +y 2 )cos<j>-2xy -\ n 
2-kJai xJ t ' L x 2 -\-y 2 — 2xycos<t> Jo 

- r * fds, 

J Aa x J t 

since inspection shows that the two values of the definite integral 
2tt and are to be used for s>r respectively, and therefore 

JAa X 

r B ' a rdr_ 
Ja'o r+z 

= r-z\og(r+z) 

L -iA'a 

This is for the radial flow from a. Adding the corresponding expression 
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for the radial flow towards b, we have finally for the complete integral 
S^-(-A% + A't + B'a-m-^ ^X^X% (4a) 
which becomes, if z = 0, 

The first line of Table I can now be filled in at once since the in- 
tegrations have shown that the first two values of k are 1 and ?; the 
next two entries are also \ since by symmetry N jr s = N jr ff = N jtr a ; 
N jc n = Njfjr = since the nadir and zenith flows are perpendicular to 
the JC flow in the horizon plane. The remaining six entries in the 
first row are for closed flows which are expressed as differences be- 
tween the flows already considered, and the corresponding k's are 
the differences of the k's for the component flows. The first column 
of k's may also be filled in, the table being symmetrical, since inter- 
changing capital and small script letters leaves N unchanged and * 
is a special case of ■*■. 

The second row of the table involves only special cases of iVjv 
and only the values \ and occur. 

From the flows included in the table nine pairs of closed flows may 
be formed having zero mutual inductances, because one of the closed 
flows of each pair has no magnetic field below the surface of the 
earth and the other closed flow includes no current above the sur- 
face of the earth, and thus there is no interlinkage of induction be- 
tween the two closed flows. The portion of Table I referred to is 
repeated in Table II, where the flows at the top are those for which 
any difference such as (>* — «) has no magnetic field below the earth's 
surface, just as {Jf—£) has no magnetic field above the earth's surface, 
as was proved above, while no flow at the side penetrates above the 
earth's surface. 
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The top row of Table II includes only values of k already found, and 
the remainder of the first column follows from symmetry and the 
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fact that a is a special case of the general flow *. Any other entry in 
Table II is now found in terms of three of the entries in this border, 
thus 

= N(jr-e)(A- a )=NjrA-Njr a -NeA+N<?a = (i-h-h)A+Ne t „ 

and we have the interesting result Ns a = 0; the remainder of the table 
follows. 

The result Nty*= — A may be readily checked directly since it 
involves only the mutual Neumann integral between straight parallel 
filaments, and by using the expanded form of the integral for equal 
filaments beginning at a common perpendicular with opposite positive 
directions 4 the result can be written down at once. 

The important difference A which is utilized in Table I may be 
expressed in the following useful forms : 

A = (-Aa+Ab+Ba-Bb) 
d 2 R 



-IS 
-If 



dSds 



dSds 

dVdv+dWdw 
R 



(5) 



= 2 ] sin \ (0 a -0i) sin [5(^1+^2) -e] ds, 

where the notation is the same as for formula (3) above. The third 
form of (5) shows that when the separation R is great the mutual 
inductance varies inversely as the first power of the separation. 

5. Mutual Inductance Between Grounded Circuits Lying on 
the Surface of the Earth 

It has now been shown that for direct currents the mutual inductance 
between grounded circuits consisting of conductors lying on the surface 
of the earth and grounded at their terminals is equal to the mutual Neu- 
mann integral between the conductors alone, since in the complete Neu- 
mann integral for the closed flows the total contribution of those parts 
which involve the ground returns is zero. For low frequencies the 
effective inductance can differ but little from the direct-current in- 
ductance, and it is therefore of practical importance to investigate 

4 Mutual Inductances of Circuits Composed of Straight Wires. Physical Review, 
5, pp. 452-458, June, 1915, formula (6). 
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the numerical magnitude of these Neumann integrals between grounded 
circuits. In order to visualize the magnitudes involved and supply 
means by which they may be readily calculated, a number of diagrams 
have been prepared for the important case of straight conductors. 5 




Fig. 4 — Mutual Neumann integral between filaments forming opposite sides of a 
rectangle with unit separation. The dashed curves show the mutual resistance 
between the filaments if the circuits are grounded through the earth and its resistivity 
is p = 2w. The numerator and denominator of the bracketed fraction on each section 
of the curve show the factors by which the vertical and horizontal scales must be 
multiplied for use on this section 



i If the two straight conductors OA, Oa start from a common point 
0, the mutual Neumann integral is shown by Fig. 2; the curves give 
the locus of terminal a for constant values of the integral when the 
other conductor OA is the unit base. The Neumann integral between 

'The necessary formulas are given in the paper loc. cit. Additional transforma- 
tions of these formulas are given in the appendix to the present paper. 
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any two straight filaments & and r having terminals A, B and a, b 
may be expressed as 



N& r = N {A a) - N iAb) - N {Ba) + N {B b) , 



(6) 



where N(Ao) stands for the Neumann integral between the two straight 
filaments, OA and Oa, beginning at 0, the point of intersection of & 
and r, extended if necessary, and ending at the terminals A and a; thus 
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Fig. 5 — Mutual Neumann integral between filaments forming opposite sides of a 

rectangle with unit separation. The dashed curves show the mutual resistance 

between the filaments if the circuits are grounded through the earth and its resistivity 

is p = 27r. This is Fig. 4, but with logarithmic scales 



the integral between any two filaments which would intersect within 
a finite distance may be readily found after reading four values from 
Fig. 2. 

In the special case of parallel filaments Fig. 2 fails; the corresponding 
curves for this case are presented by Fig. 3, which assumes a unit base 
filament and a parallel filament starting at a point on the left hand 
perpendicular to the base. Differences will give the general case 
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of parallel filaments which do not start at a common perpendicular, 
which may thus be derived from Fig. 3. 

The mutual Neumann integral between any two parallel filaments 
may also be obtained by means of the formula 

Ns?r = M - N {Aa) + N {Ab) + N (Ba) - N iBb) ], (7) 

where N {Aa ) now stands for the Neumann integral between the pro- 
jections oi Aa on SP and >-, extended if necessary, the projections 



-■Y. , -°* -¥: , ,-?,* °, 3 , -V. -Of ■ ■ .9 , P.'- , ?? "l 3 °, a °f V "ft I . "fl . . P,\ I .'j* . 




-10 -Q9 -Ofl -07 -08 -0.8 -0.4 -0.3 -02 -0.1 



0.1 02 03 0/> 0.5 0.0 07 03 OB 10. 



Fig. 6 — Contour lines of the mutual Neumann integral between a counter-clockwise 
small loop on the surface of the earth, per unit area, and a straight grounded filament 

AB of unit length 

having the same or opposite positive directions in agreement with 
.2? and r. This formula for the mutual Neumann integral presents the 
advantage of requiring only a single entry diagram, which is supplied 
by Fig. 4 and on a logarithmic scale by Fig. 5. 

The mutual inductance may be required between a small, closed 
loop lying upon, but insulated from, the surface of the earth and a 
straight grounded conductor. The value depends upon the location, 
area and assumed positive direction around the loop, but is inde- 
pendent of the shape of the loop. Contour curves for the mutual 
inductance per unit area of the loop are given by Figs. 6 and 7; the 
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positive direction around the loop is counter-clockwise; the straight 
grounded conductor AB of Fig. 6 is of unit length while in Fig. 7 
grounded terminal B alone appears, terminal A being at an infinite 




Fig. 7 — Contour lines of the mutual Neumann integral between a counter-clockwise 
small loop on the surface of the earth, per unit area, and a straight grounded filament 

of infinite length 

distance to the left. These curves show the vertical component of 
the magnetic field due to unit current in AB. The formulas employed 
for AB = \ and infinity, respectively, are B 



d 2 N 
dxdy 

d 2 N 
dxdy 



_ 2yf>i+r 2 ) = ri+r 2 Jl 
rir 2 [(r, + r 2 ) 2 -l] r x r 2 \ (, 



-(r,-r 2 ) 2 



ri + r,)*-l' 



—. r-Z- — r = - tan \0 2 . 
riKXi+r-i) r 2 



(8) 
(9) 



Large loop mutual inductances may be calculated either by inte- 
grating the value of —d 2 N/dxdy over the loop or by integrating the 
value of dN/ds around the boundary. If the boundary may be ap- 
proximated to by a broken straight line, the curves of Figs. 2 and 3 
may be employed. 

6 These formulas may be derived by differentiating (8) of the paper loc. cit., 
dN/dx = 2 tanh - l[4A/(Ai+fB)], with respect to y. 
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6. Mutual Impedances for Conductors Lying on the Surface 

of the Earth 

In order to arrive as directly as possible at a concrete numerical 
idea of the magnitudes and angles occurring in the mutual impedances 
encountered in engineering work, we may advantageously start with 
the following specific constants: 

Base Length {AB, OA or A a) 1 Mile, 

Frequency of the Alternating Current 1 Kilocycle, 

Resistivity of the Earth per Centimeter Cube . 1 Megohm, 

which are of the right order of magnitude and make the factors 

p/(2tt AB) =10 6 /(2tt 0.1609 X10 6 ) =0.989, 
27t/X^J3X10- 9 = 2t 0.1609 = 1.011, 

which are both equal to unity within about 1 per cent., so that the 
approximate resistance and reactance components of mutual im- 
pedances may be read directly from Figs. 1-8 without applying multi- 
pliers. On the other hand, when mutual impedances are required 
for other lengths, frequencies and specific resistances, the correcting 
factors are readily applied. The tangent of the angle of the mutual 
impedance is proportional to the frequency, to the square of the 
linear dimensions of the circuits and to the reciprocal of the earth's 
resistivity. 

Grounded circuits separated by a distance large compared with 
the dimensions of the circuits have a mutual impedance with negligible 
resistance component since ultimately this component varies in- 
versely as the cube of the distance by (3), whereas ultimately the 
mutual reactance varies inversely with the distance. 

For two parallel grounded conductors separated by one mile and 
forming opposite sides of a rectangle, the two components of the 
mutual impedance are shown by Fig. 5 for the assumed constants 
(approximately a kilocycle and a megohm). The resistance com- 
ponent of the mutual impedance is then always less than the reactance 
component; when the rectangle becomes a square, the mutual im- 
pedance angle is tan -1 1.595 = 57.9°. Reducing the frequency to 
0.627 kilocycles or reducing the side of the square to 0.792 miles or 
increasing the resistivity to 1.595 megohms would reduce this angle 
to 45°. 

Consider two straight grounded conductors AB and ab, the latter 
distance being small compared with the other dimensions of the 
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system and assume that while ground a is fixed, ground b is rotated 
about a in a circle of fixed radius ab. This will vary the mutual 
impedance between the two grounded conductors. The reactance 
component will vary as the cosine of the angle between ab and AB. 
The resistance component will also vary as the cosine of an angle, 
measured, in general, from a direction other than AB. The maxi- 
mum resistance component will also differ, in general, from the maxi- 
mum reactance component. The locus of the mutual impedance 
obtained for all positions of ab will be an ellipse. The ellipse becomes 
a straight line when ground a lies on the bisecting normal of AB, 
for then the direction giving the maximum resistance component is 
parallel to AB and thus the same as for the maximum reactance com- 
ponent. The straight line limit is also obtained when ground a lies 
on the prolongation of AB in either direction, for the resistance 
component then has its maximum value in the direction opposite to 
AB. The maximum resistance component will be perpendicular 
to AB at the points on Fig. 1 where the C/I contour, if drawn, would 
be vertical. The locus of these points is 

y 2 = (x 2 -b*)VH{x+b)W+(x-b) 2 / 3 ]. (10) 

At remote points on this locus y = ± V2 x and the maximum resist- 
ance is negligible compared with the maximum reactance since the 
circuits are widely separated, while in the neighborhood of A and B 
it is the maximum reactance which is negligible compared with the 
maximum resistance. At some intermediate point the two maxima 
are equal, and, since they occur for directions differing by 90°, the 
elliptical impedance locus becomes a circle; and the mutual impedance 
between the two grounded circuits does not change in magnitude 
as ground b is rotated about ground a. For the assumed constants 
(1 mile, 1 kilocycle and 1 megohm) this point lies at distances of 
1.562 and 0.939 miles from the two terminals A and B, and its four 
possible locations are shown by the four small crosses on Fig. 1. 

If a is rotated counter-clockwise about AB, the direction giving the 
maximum resistance component also rotates counter-clockwise, making 
two complete revolutions, while the ground a makes one revolution 
about AB. 

7. Equivalent Ground Plane 

To a first approximation the direct-current mutual inductance 
between two straight conductors AB and ab, forming opposite sides 
of a rectangle on the surface of the earth, at a separation Aa which 
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is small compared with the length AB is, by (25) of the appendix, 
neglecting the first and higher powers of 1/s, 7 

*AB 

N = 2AB\og^- i —=2AB\og ™ AB . (11) 

b Aa & Aa ' 

This expression has the form (2 I log s/r) of the commonly employed 
mutual inductance formula for two long parallel conductors, each of 
length /, separated by distance r, the common return being a perfectly 
conducting earth in which the image of each conductor is at the 
distance s from the other physical conductor. For our direct-current 
case, therefore, the effective distance to the images is about j of 
the length of either grounded conductor. Since this distance is by 
assumption large compared with the distance between the conductors, 
the images are approximately at this same depth below the actual 
surface of the earth, and the hypothetical perfectly conducting earth 
would be at one-half this depth, or \AB. The effective image dis- 
tance is necessarily directly proportional to the dimensions of the 
grounded circuits and independent of the earth's resistivity because 
the shape and relative distribution of the lines of flow are independent 
of the resistivity and of the length of the grounded circuits. Inspec- 
tion of Fig. 1 shows that somewhat over \ of the return flow attains a 
distance \AB from AB, while the remainder of the current remains 
closer to the grounded conductor. 

It may be inquired what would be the effect of confining both return 
currents to a thin uniform conducting layer on the earth's surface, 
so that they become horizon return flows. For the closed flows 
(Jf—je) and (■«—■*) in general and the particular flows (&-■&) and 
(r— a), where ^P and r are close parallel straight conductors, 

N(jr-j?){jt-/>) = Njr*-NjrA-Nj? J ,+NjeA 

= Njf J r-A; 
N(SP-^)( r -A) = Nsp r -2Ab + 2Aa 

*AB 
= 2AB log t—. =2AB log ^j^. (12) 

7 If the term \/s = Aa/AB of the expansion is retained the equivalent ground plane 
has the depth (AB+Aa)/e and thus becomes deeper as ab is moved away from AB. 
But the equivalent ground plane may be kept fixed at the distance AB/e from AB 
provided it is tipped at the angle sin - 1 2/e=47° so that ab moves away from the 
ground plane as it moves away from AB. If it were worth while, still closer ap- 
proximations might be secured by using a perfectly conducting cylindrical earth of 
suitable cross-section. 
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Thus, the assumption that the return current is confined to the earth's 
surface does not change the order of magnitude of the effective image 
distance, but reduces it from about f to about \ of the length of the 
exposure. For space returns the effective image distance is 

2AB/e zl ' i = 0Ai§AB. 

Now take another practical case by assuming that the conductor 
ab is of negligible length compared with its separation r from the 
parallel conductor AB and the separation is, in turn, negligible com- 
pared with the length AB. The formula for dN/dx given in footnote 
6 shows that the required inductance depends only on the ratio 
AB/iri+rz) and is thus constant upon an ellipse. Equivalent ex- 
pressions in logarithmic form are 

N = 2 ab log (2 cos §0, sin fo^Ip) (13) 

= 2aMog (^jM, (14) 

or approximately 

AB 
N = 2 ab log , if ab is opposite the midpoint of AB, (15) 

N=- 2 ab log — *— , if ab is at distance r 2 beyond B. s (16) 
2 L ti — I 

Thus, from (13) the effective image distance is never greater than 
twice the geometrical mean distance from ab to A and B. Its maxi- 
mum value is approximately AB and occurs when ab is opposite 
the midpoint oi AB. Its minimum value is approximately Vr 2 iB 
and occurs when ab is at the distance r% from A or B in the prolonga- 
tion of AB. This makes the inductance one-half of what it is at the 
symmetrical position. Thus, wherever ab is placed, its mutual in- 
ductance with AB lies somewhere between 50 per cent, and 100 per 
cent, of the mutual inductance, due to an effective image distance 
AB, ab remaining always parallel to AB and the locus of ab being a 
rectangle with semicircular ends of radius r 2 and centers A and B. 

8. Mutual Impedances-of Grounded Circuits Which Depart 
from the Surface of the Earth 

Consider a system of conductors following any paths in space and 
insulated from the earth except at two grounding points on the sur- 

8 The shortest distance from ab to AB might have been designated by a single 
letter in place of using y, r and r 2 in formulas (13), (15) and (16). , /. ^ 
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face of the earth. Any flow of current through this system of con- 
ductors may be divided into elementary filaments each of which is 
made up of segments beginning and ending at the earth's surface 
and not crossing the earth's surface between these terminals. Ground 
each segment at both ends. Let & and 1/ designate segments having 
terminals A and B, W (for example, an open wire circuit) never 
going below the surface of the earth and W (for example, an under- 
ground cable circuit) never going above the surface of the earth. Add 
the underground flow W to Table II at the foot of the left-hand 
column and, proceeding as before, 

= N{9f~ jr)(.r-«) = N&M-Ntta'+ib = N&,+N2/*-h*t 
= N{2/- X) (A - a) = Nt/a - Ni/ a - ?A = Nt/a + Nffi, - f A , 
or NT/a - JA- Afar = f A - Nt/a, 
and similarly for the flow W above ground, 

Nw e = \A + Nwn = Nwj, - |A. 
Hence the three cases which may occur give 

N(W-£){w-r) = N7i/u,-N?i/„ - Ntyu, 

= Nwu,- Nwa- Nj?u,+2A, J 



(17) 



N{9/-£){ u - tf )=N?/u+N9/*+N3ru 1 

= N?/ u +N2sa + Nj? u -2&, J 

N{W-iS){u-<,) = N?£/u-Nw„-\-Njru 1 

= Nw u - Nwa + Njffu. J 

The importance of these equations lies in the fact that the earth 
return flows are replaced by the simpler nadir, zenith and horizon 
return flows. // the conductors comprise only broken straight filaments, 
making any angles with each other and the earth, the required Neumann 
integrals, if we use the expressions involving the nadir and zenith returns, 
are the known expressions between straight filaments. If the conductors 
lie in horizontal planes with vertical ground connections, it is con- 
venient to employ the expressions involving the horizon return flows, 
since the required integral is (4a) derived above. Formulas (33)-(41) 
of the appendix are the resulting formulas for the three general cases 
and for a number of important special cases. 

In general we may say that the effect of changing the height of one 
or both conductors by an amount which is small compared with the 
length of the conductors will be relatively small, since the effective 
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image distance has been shown above to be of the order of thejength 
of the conductors. To illustrate this fact, in Fig. 8 four dotted curves 
have been added to the curve of Fig. 5 showing the mutual induct- 
ances of the two parallel grounded conductors when they are in the 



100 




Fig. 8 — Mutual Neumann integral between grounded filaments forming opposite 

sides of a rectangle with unit separation, on the surface of the earth (from Fig. 5) 

and between the same filaments when one or both of the filaments is raised above 

or depressed below the surface of the earth as shown by the insert 



four positions indicated by N\, iV 2 , Ns and N* on the insert of Fig. 8, 
calculated by formulas (28)-(30) of the appendix, which are special 
cases of (35), (36), (39) and (40). When the conductors are long, the 
relative change in the mutual inductance is small. Depressing the 
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Fig. 9 — Comparison of measured and calculated mutual impedances between the 

indicated circuit pairs. Each arrcw extends from the measured impedance located 

at its tail to the calculated impedance at its head. The location of and positive 

directions for these circuits are shown by Fig. 10 



24 



BELL SYSTEM TECHNICAL JOURNAL 



wires reduces the inductance; the curves show that in the case of 
JV 4 the inductance passes through zero and is reversed in sign when 
5 = 1.560. 

When the departure of the circuits from the earth's surface may be 
neglected, all terms, but the first, on the right-hand side of formulas 
(17)— (19) drop out, and each reduces to the simple, fundamental 
grounded circuit formula (4). 

9. Comparison of Theoretical Results with Measurements 
at 25 and 60 Cycles 

Fig. 9 shows, by means of arrows, the impedances which must be 
added to each of a large group of measured 25-cycle mutual im- 
pedances to obtain the results calculated by means of the preceding 
formulas, on the assumption that the earth has a uniform resistivity 
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Fig. 1C — Location of test circuits, the arrow-heads showing the positive directions 
in the conductors. When circuits are combined in series, with the removal of 
intermediate grounds, the new circuit designation is shown by the equations. The 
test circuits A, m and t are large metallic loops from which all grounds have been 
removed. The horizontal and vertical displacements of the conductor by a few feet, 
which render the indicated equalities for -^and n only approximate, were allowed 
for in determining the calculated results for Fig. 9. The grounds of the capital 
letter circuits were isolated sections of single track about one mile in length; the 
midpoint of the section is shown as the effective ground but it may have actually 
been displaced and have varied with the moisture of the road-bed 
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of 0.5 megohm per centimeter cube. The measurements were not 
made for the purpose of this comparison, for which they are not well 
adapted, but they do give both components of the mutual impedance, 
which is the absolutely essential requirement. The geometrical 
irregularity of these circuits is shown by Fig. 10. This was com- 
pletely allowed for by making detailed computations after substitut- 
ing an approximately equivalent broken line for each circuit. The 
variability in the earth's resistivity with location, depth and changing 
moisture content on different days could not be allowed for. The 
effect of buried gas and water pipes and of other grounded conductors 
was also necessarily neglected. 

The direct-current theory leaves but one arbitrary constant at our 
disposal after the frequency and the geometry of the circuits have 
been fixed. This constant is the earth's resistivity. By trial it 
was found that 0.5 megohm gave a good average agreement between 
the calculated resistance components and the entire set of measure- 
ments, only a part of which is included in Fig. 9. The individual 
discrepancies are large but are not so large as to be disconcerting, 
considering the variations in effective earth resistivity from place to 
place and from day to day during the progress of the tests. 

The calculated reactance component of the mutual impedance 
based on the direct-current mutual inductance is independent of the 
earth's resistivity and is uniquely determined by the frequency and 
geometrical relations. Even a general agreement between the cal- 
culated and the measured reactances is significant and Fig. 9 shows not 
only this, but also a great many good individual agreements. The 
outstanding discrepancies for circuits C and <? are systematic, and are 
apparently to be explained by the effective grounding of these cir- 
cuits at some other points than the midpoints of the track sections. 
On the basis of this comparison, it appears that the direct-current 
theory proves itself adequate to give an approximation to the actual 
mutual reactances, provided the linear scale and the frequency in- 
volved do not greatly exceed those of these tests. 

Measurements were also made at 60 cycles. The resistance com- 
ponent remained roughly the same as for 25 cycles; the reactance 
component doubled as shown by Fig. 11; each component therefore 
agreeing approximately with the results which would obtain if the 
direct-current distribution is maintained. 

Other comparisons have been made with the same conclusion, but 
tests should be made, throughout a range of frequencies, at a locality 
where it is known that the conductivity of the upper layer of the 
earth's surface is reasonably uniform so that the effect of the lower 
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layers may be determined. Small scale models, having the proper 
propagation constant, could advantageously be used in determining 
the alternating-current impedances for uniform earth resistivity or 
any assigned distribution of resistivity. 
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Fig. 11 — Comparison of measured mutual reactances at 25 and 60 cycles. A right 

line with the slope 2.0 fits the point somewhat better than a line with the slope 2.4 

which corresponds to the ratio of the frequencies 



10. Summary 

Formulas for the direct-current mutual resistances and mutual 
inductances for grounded circuits on the assumption of uniform 
earth resistivity have been derived and useful diagrams prepared. 
The applicability of these results as a first approximation to many 
practical alternating-current cases has been shown. 

If, as I hope, this paper is free from ambiguities and errors, it is 
due to a thorough revision by Mr. R. M. Foster; and I am indebted 
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to Miss Frances Thorndike for the accuracy attained in the numerous 
curves of Figs. 1-5 and 8, which should make them of practical value 
in numerical calculation. 

Mathematical Appendix 

Additional mathematical results which have been employed in 
connection with the figures and discussion of this paper are brought 
together below for convenient reference. 

Formulas for Fig. 2 

2s = P +d+l,tiOA = l, 0a = P , 

d 2 = l+p 2 -2pcos0, 
sin0 = dsin (0+0), 
sin = p sin (0+0), 

COS =(25-1)- 25(5 -l)/p, 



j8^ 

/ / 


a 





\ ' "A 

Fig. 12 



A3 



=f(p,e)= P f(p-\d) 



= cos (log -i- +p log -^ ) (20) 

op o X 



= COS V 



. tan £(0+0) . . 1 ,, m 

l0g tan TO +Pl ° g tan£0tani0 > (21) 



= 2cos0log^=2cos01og(l + l/sin£0),ifp=l, 

= £plog|^+|pMog2±- p ,ifd = l, 
2 — p p 

= cos0|^plog-^ -(l-p)log(l-p)J 

= cos/3[plogp-(l+p)log(l+p)]- i ^p-^^ 

p(ll + 19 P +llp 2 ) 
" r " 96(1 +p) 3 p T ' " ' * 

= _ 2 _ d ,-[tf+(i-p)iog(i-p)j/2P_|_ ^ (22) 

Vp 
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21og(l+p)-p/(.l+p) 



when0 = 7r) (23) 



= 0, 1.564, co at p = 0, 1, co, which checks the 
sharp curvature of the curves which cross 
the axis just to the left of the origin. 

Formulas for Figs. 4 and 5 
s=AB/Aa, 



NjPr 



(25) 



^ = 2[slogO> + vV + l) + l-\A 2 -r-l] (24) 

= 2,[log2,-l + ^-l +3 4 ? -i + ...] 
= 25 log [ 2 i(l+i+2_- I l- 3 -J- T+ .. .)] 

-'E 1 a+a-w+ ■"••]• (26) 

^^)(y)= 2 -vfc' / ' "' (27) 

^J = ^T+2 1og(5 2 +l), (28) 

iV 2 (orAM = i^- ±2[1 x (1 + V 7q^- ) _V,^+l+^+l]. (29) 

yla ,4a 

N,_N&r. 2 f lo (5 2 +l) (V^+i) 4 
4a Aa L° g (i + vV+2) 4 

+4(V^+2-V^+T-Vl+l)]. (30) 

Formulas for the Mutual Inductance Between Any Flows in Two Hori- 
zontal Planes Grounded by Vertical Filaments 

Let the arbitrary flows be <X" and *' between points A', B' and a', b' 
in the two horizontal planes, grounded by vertical filaments connect- 
ing these four terminals with the points A, B, a, b on the surface of 
the earth. In order to indicate briefly which of these eight points 
are involved in each term of the result, we imagine a vertical line 
which cuts the horizontal planes in the points P', p', P, p, where P 
and p are the same point on the surface of the earth, since the non- 
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primed points are all in this plane, and we agree that A' or A occurs 
in a term, according as P' or P is found in the subscript of the symbol 
A or T used to designate the term, where 



A P - p -=-A'a'+A'b'+B'a'-B'b', (31) 

{A'b'+P'p'){B'a'+P'p') 
Tpp ' ~ IOg {A'a'+P'p'){B'b'+P'p') ■ {6Z) 



In these expressions every distance between points, such as A'b' , 
P'p', is a positive quantity. The formulas below are perfectly general, 
but require the assignment of the capital letters -X"' , A', B', P' to the 
upper plane when both flows are above the earth and to the lower 
plane when both flows are below the earth. They are most readily 
checked by employing formulas (4a) and (24) in formulas (17), (18) 
and (19). The results show that the mutual inductance is equal 
to the Neumann integral between -**' and *■' augmented by terms 
which depend only upon the arithmetical distances between the eight 
points A', B', a', b', A, B, a, b. 

N{w-f)(u,-e) = Njir>*.'+P'p'rp- P '+2Pp'rpp-Ap V +App, 

where P'p>Pp', (33) 

= ^ V +2Z logoff], if i* and,- are 

both at height Z, (34} 

= iVjrV + 2Z log (1+5 2 ), if A'B'b'a' is a 

horizontal rectangle and AB = s(Aa), (35) 

= A^rV + 2Z[logi(l + VT+?) + 1-Vr+^ + s], 
if A'B'b'a' is a vertical rectangle with one 
side on the earth and AB = s(A'a)=sZ, (36) 

Nw-£ ) (u-*) = NJC'J + P'P'Yp-p' ~ 2P'pr P ' P - 2Pp'(r Pp - - T PP ) 
-A Pp - + 2Ap- p + 2Ap P '-'SAp p , 
where P'p>Pp', (37) 

-isa * -£ log ( ^ 6)(5a)( ^ a ' +Z )-->(56'4-Z) 2 
+-i(Ap> p -Ap p ), if P' and p' are both 
at distance Z below the earth, (38) 
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(/+Vi+* 2 +/ 2 )< 

+8 (Aa) (Vl + * 2 + * 2 + 1 - Vl + s 2 
- Vl+t 2 ), xlA'B'b'a' is a horizontal 
rectangle at the distance Z = t(Aa) 
below the earth and A B = s(A a) , (39) 

=iVjrv-2Z[iogKi+Vr+^)-\/i+^+i+^], 

if A'B'b'a' is a vertical rectangle with 
one edge at the surface of the earth and 
AB = ab^s(A'a)=sZ, (40) 

N{W-£)(u-.) = Njr' J r'+P f p'r P >p'-2Pp , rpp—Ap p >+2Ap l> >-Ap p . (41) 



